Water waves over a time-dependent bottom: Exact description for 2D potential flows 
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Two-dimensional potential flows of an ideal fluid with a free surface are considered in situations 
when shape of the bottom depends on time due to external reasons. Exact nonlinear equations 
describing surface waves in terms of the so called conformal variables are derived for an arbitrary 
time-evolving bottom parameterized by an analytical function. An efficient numerical method for 
the obtained equations is suggested. 
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I. INTRODUCTION 



The theory of water waves now is among the most ac- 
tively developing branches of the hydrodynamic science 
i^e^or instance, Illllllllli,|i|l3,|lll[ll 
ll^ Im . Hsl IT^ . and many references therein). In par- 
ticular, two-dimensional (2D) flows of an ideal incom- 
pressible fluid with a free surface are very attractive for 
theoretical study, since they possess many essential fea- 
tures in common with real water waves, and simultane- 
ously are relatively simple for analytical and numerical 
treatment. A powerful mathematical support for the 2D 
theory is provided by analytical functions and the corre- 
sponding conformal mappings. The original idea of us- 
ing conformal mappings was intended for waves on the 
deep water or over a static horizontal bottom: a ge- 
ometrically simple region occumed by restin g flu id (it 
is the lower half-plane H S [li [H El Ellll, or a 
horizontal stripe 12J) is mapped onto the flow region 
with a disturbed time-dependent free surface, so that 
the real axis is transformed into the moving boundary. 
With such conformal "straightening" of the free bound- 
ary, it is possible to obtain an exact analytical description 
for surface waves. The corresponding (l-l-l)-dimensional 
nonlinear equations of motion are very convenient for 
numerical simulations by spectral methods. Important 
achievements in this direction were reported in works 
H E in III 111 111 m, including analytical and numer- 
ical results. A generalization has been recently made to 
the case when the bottom shape is still static but strongly 
inhomogencous in space |16|. 

The purpose of present paper is to derive exact equa- 
tions of motion for the free boundary in 2D ideal 
flows, valid with arbitrary non-stationary bottom profile 
(Sec. II), and to suggest a numerical method for solving 
the obtained highly nonlinear equations (Sec. III). This 
problem, besides being interesting by itself, is also impor- 
tant for studying such phenomena as tsunami generation 
by earthquakes and other similar processes in natural and 
laboratory conditions. 



II. EXACT NONLINEAR EQUATIONS 

Let us start our consideration with basic definitions. 
We will describe a nonstationary flow region 7)^ in (x, y)- 
plane (with y-axis up-directed) by an analytical function 
X + iy — z(w,t), where t is the time, and the complex 
variable w = u + iv occupies in (u, w)-plane the flxed 
inflnite horizontal stripe : —oo < u < +oo, < v < 
1. Position of the free surface dvi'^'' will be given in 
parametric form by the formula 

X'^''\u,t) + iY(''\u,t) = Z(^)(u,i) = z{u + i,t), (1) 

while the bottom profile 9Pi''-'will be determined by 

X^''\u,t) + tY'-''\u,t) = Z^''\u,t) = z{u,t). (2) 

Inasmuch as the velocity potential Lp[x,y,t) satisfies the 
Laplace equation (f^x + '^yy = 0, in "Dz the analytical 
complex potential 4'{z,t) — Lp{x,y,t) -\- iO{x,y,t) is de- 
fined, with 9{x,y) being the harmonically conjugate for 
^p{x, y). We deal below with the complex analytical func- 
tion (f){w,t) — (j){z{w,t),t) defined in I?^ and taking the 
boundary values as written below: 

0(it + i, <) = $(")(«, i), 0(u,t) = $('')(w,t). (3) 

The velocity components are determined by the usual 
relations 

Vx - iVy ^ (fix - i^y = d4>/dz = 4)' {w,t)/ z' {w,t), (4) 

where (. . .)' denotes derivative over the complex variable 
w. 

Now we proceed to equations of motion. First of all, 
it is the dynamic Bernoulli equation on the free surface. 



Re(^t(z,t))-h^^^^^+5lmz = 0, 



zedVi^\ (5) 



where g is the gravitational acceleration. This equation 
can be easily re-written in terms of the conformal vari- 
ables as follows: 
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4, ) + + g Im Z(^) = 0. (6) 
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Two other equations of motion are the kinematic con- 
ditions on the free surface and on the bottom, expressing 
the fact that the boundary motion is determined by the 
normal component of the velocity field. In the conformal 
variables these conditions take form 



Imfzi'^zW^l =-Im$?) 



(7) 



(8) 



where Z denotes complex conjugate value. 

We should take into account that ^'^"^ {u, t) and 
^^^^{u,t) are not mutually independent, but they are re- 
lated to each other by a linear transform. Indeed, the 
following formula for analytical continuation is valid: 



(f>{w,t) 



(b) ikw 



dk 

2^' 



(9) 



Here the Fourier transform <I>*^''^ (m) ^ ^j. is defined in 
the usual way: 



(10) 



Therefore at the free surface, where w = u + j, we have 



dk 



27r 



(11) 



It is a proper point here to introduce some linear op- 
erators, necessary for further exposition. These linear 
operators S*, R, and T = are diagonal in Fourier 
representation: 

Sk — l/cosh(fc), Rk = itanh(fc), Tk = — icoth(A;). 

(12) 

It will be convenient for us to write {u, t) in the form 
(u, t) = Sij{u, t) - i{l - iR)dZ^f{u, t), (13) 

with ip{u,t) and f{u,t) being some unknown real func- 
tions. Then for <^^^\u^t) we will have the formula 

$(")(«, = [l + iR)Hu,t) - iSdZ^ f{u,t). (14) 

Later we use in equations the complex function 

*(w,t) = (l + ii?)i/.(u,t). (15) 

Now it is necessary to look at the conformal mapping 
z{w,t) with more attention. Since the bottom motion is 
assumed to be prescribed, this function has the following 
structure (compare with [16|): 



z{w,t) = Z{C{w,t),t), 



(16) 



where a known analytical function Z(C, i) determines a 
conformal mapping for the upper half-plane of an in- 
termediate complex variable C onto the infinite region 



above the bottom in z-plane (or, at least, the mapping 
Z((^, t) should have no singularities within a sufficiently 
wide horizontal stripe above the real axis in ^-plane, so 
there should be enough "free space" for large-amplitude 
waves in z-plane). The intermediate analytical function 
C{w,t) takes real values at the real axis. Therefore we 
may write 



cosh(fc) 27r' 



(17) 



where ak{t) is the Fourier transform of a real function 
a(u,t). On the bottom C,{u,t) — Sa{u,t), thus the para- 
metrically given curve 



z'-''\u,t) = ZiSa{u,t),t), 



-oo < U < +00, 



(18) 



coincides with the prescribed curve X + iY = 
Z{s, t), —oo < s < +O0 (these two curves differ from 
each other just by parameterization). At the free surface 

C(u + i,t)= ^(u, t) = {l + iR)a{u, t), (19) 

Z'^''\u,t) ^ Z{C{u,t),t). (20) 

Since our purpose is to derive equations of motion for 
the unknown functions 5'(it,t), f{u,t), and ^{u,t), we 
have to substitute expressions ((T^ . ifTH) . lfTS|l . ((T^ . and 
(1201 into Eqs.©, 0, ©■ Eq.® is simple and does not 
need any trick. Eq.© (divided by IZ^*-*!^) takes the form 



Im ( 1^ 



Im 



-Im ^If, 



Sf 



(21) 



From here we express (see Eg . (1231) below and compare 
with ^^) and substitute proper expressions into Eq.®, 
in order to fi nd i )f Simple transformations, similar to 
those in work [13, lead to the following equations: 



/= \TL\[Zt{s,t)Zs{s,t)Sau 



s—Sa 



6 = (T+i) 

{f+l 

+ {l + iR) 



"lm*„ - Sf 

*«^t(?,t) 



™U?(c,i)en 



Ret (22) 
,(23) 



Im 



Re 



zdiMu 
'2\z^{iMu? 



-<7lmZ(C,i) 



(24) 



This system of equations is the main result of present 
paper. It provides exact description for potential water 
waves over a time-dependent bottom determined by an 
analytical function Z{Q,t). It should be noted that cap- 
illarity effects can be easily included here by adding the 
term 

a\d,,Z{au))\-^\u,[{dufZ{^{iL))/d^Z{i{u))] 
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FIG. 1: (i) Initial bottom profile and horizontal free surface. 
The velocity field is everywhere zero. 



(proportional to the surface curvature) to the gravity 
term —gln\Z(^). 

The above equations are valid not only for infinitely 
extended configurations (when the fluid covers the whole 
bottom), but also for the flows with one or two non- 
stationary shore edges. However, asymptotic properties 
of the intermediate mapping C(w^j i) £^re different in these 
cases. So, with two edges C(u,t) C±{t) as u ^ ±oo, 
while for edge-less configurations the corresponding lim- 
its are Ci^jt) ±oo. 



III. NUMERICAL METHOD 

Equations though look quite compli- 

cated, in many cases are convenient for numerical simu- 
lations, as described below. 

Instead of the function ^{u,t), we deal with another 
function (similar to that proposed in [T^ITH ). 



f3{u,t) = ^u{u,t)/^u{u,t), 



(25) 



because this choice results in better numerical stability 
(concerning the problem of optimal variables see 0,^3)- 
Let us consider a class of space-periodic flows (with- 
out shore edges!) over a;-periodic profiles of the bottom. 
With a proper choice for the length and time scales, we 
may write g = 1, and Z{( + 27r, i) = 27r -I- Z{(, t). As it 
was pointed in |la | , in space-periodic flows the variable u 
comes into solutions in the combination -d — ua{t), where 
a real function a{t) depends on time because the operator 
T is singular on small wave- numbers. Besides a{t), solu- 
tions are determined by two other real functions, p('(5,t) 
and x(t^, ^), both having period 27r on variable: 



H-oo 



p{i!),t) = ^ p„(t)exp(iTOz?), p_„(t) = p„i(i),(26) 
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FIG. 2: (i) t = 2: The cavity on the bottom has been almost 
formed. 



-t-oo 

X(i?,i)= Xmit)exip{imd), X-mW = XmW,(27) 



m— — oo 



-|-oc 



2/5m(i) exp(im??) 



r?i— — oo 



+ exp(2ma(t)) 
d + ia{t) + {l + iRo,)p{'0,t), (28) 



+ 00 



2 Xm(^) exp(^m^?) 
-t- exp(2mQ;(i)) 



(I + »R„)x(i9,t). 



(29) 

The linear operator Rq is diagonal in the discrete Fourier 
representation: Rq(to) = itanh(Q;m). Let us introduce 
the following quantities: 



F= Im 



C/=(T„ + ») 



Im (/3e')-S„F 



Im 



(30) 
(31) 



M.^S^+InrZ(,,.)-Reff##V32) 



2|^c(e,t)ep 
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FIG. 3: (i) t = 10: The right-propagating group of waves is 
approaching the shallow region. 

where Sa{m) = 1/ cosh(ar7i), ^' = d^^ = l + {l+iRa)d^p. 
The hnear operator Tq, is regular. In the discrete Fourier 
representation it is defined as follows: 

Ta{'m) — — icoth(am), m 7^ 0; 

= 0, m = 0. (33) 

Equations of motion for the real functions a{t), p{'d,t), 
and xi"^, t) follow from and have the form 

a{t) = Im([/)di^, (34) 

p{^,t) = -ReiUa, (35) 
X{^,t) = -Re(^{//3' + l(l + iR„)a^I^^ , (36) 

where /?' = d^f} = (1 + if^a)d^X- 

The above system (|34|) - (I36|I can be efficiently simulated 
by the spectral method, since the multiplications are sim- 
ple in representation, while the linear operators Sq, 
Rq, and Tq (also i9-differentiation) are simple in Fourier 
representation. The integration scheme can be based, 
for instance, on the Runge-Kutta 4-th order accuracy 



FIG. 4: (i) Breaking leading wave at t ^ 13.18 

method, similarly to works [Tsl Il6j| . Efhcient subroutine 
libraries for the fast Fourier transform are now available. 
In the numerical experiments reported here, the FFTW 
library was used JJj. The employed Fourier harmonics 
had numbers m in the limits —16000 < m < 16000 (the 
size of arrays for the Fourier transform was even larger, 
N — 2-^^), thus the obtained solutions are very accurate 
(compare with 16]). 

A. First numerical experiment 

In the first numerical experiment [referred to as (i)], 
the shape of the bottom was determined by analytical 
function Z{(^,t) = S'*-'(C — iao,t), where ao = tt/IOO w 
0.0314, and function B^^\q,t) is expressed as follows: 

/ N . , I i sin Q + \/e + cos^ q \ 
B^'\q,t) - q + iAlnl 

+ -^exp[-C7(l + cos(<z-J))]. (37) 

Here the parameters are A = 0.6, e = 0.01, A = —0.008, 
r = 1.0, C = 500.0, and 5 = 0.167r. The time- 
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FIG. 5: (ii) t = 5 



FIG. 6: (ii) t = 10 



independent terms in the above formula describe the bot- 
tom profile as shown in Fig^ with deep and shallow re- 
gions. The time-dependent term in B{q, t) corresponds to 
formation of a cavity on the bottom. The initial values 
were taken jo(t?,0) = 0, x(i?,0) = 0, and a(0) = ao, thus 
giving the horizontal free surface and zero velocity field 
at t — 0. The corresponding numerical results are pre- 
sented in Figs lll4l in terms of the dimensionless quantities 
z = x + iy and V = [(3C - iSaF]/[Z^{^,t)e] = - iVy. 
This numerical experiment can be viewed as a rough 
model for a tsunami wave generation, subsequent prop- 
agation, and final breaking. Indeed, the inhomogeneous 
bottom displacement is an analog of an earthquake re- 
sulting in two opposite-propagating groups of surface 
waves. When one of the groups first approaches a shallow 
region, the crest of its leading wave becomes higher and 
sharper, and finally singular. In real world instead true 
singularity one can observe vortices, splash, and foam on 
the crest. Further treatment of this situation is not pos- 
sible within present ideal 2D theory, since 3D effects and 
dissipative processes become important. 



B. Second numerical experiment 
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In the second numerical experiment [referred to as (ii)], 
the bottom profile was determined by function Z{(^,t) = 



FIG. 7; (ii) t = 10.117 
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FIG. 8: (iii) t = 5 
S("nC-iao,t), where 



FIG. 9: (iii) t = 10 



sin q + y^e + cos^ 



VT+~e J ' 

(38) 

with ao = 7r/100, e = 0.01, A = 0.1, A = 0.25. The initial 
values were 0) — 0, x('?, 0) = 0, and a(0) — ao, thus 
corresponding to the horizontal free surface y = 0, the 
horizontal bottom y — — ao, and zero velocity field at 
t = 0. At later times, some pieces of the bottom rose 
and some pieces went down. Results of this simulation 
are presented in Figs l5l7l In this case, a singularity forms 
in the places where the fluid flows from the shallow region 
to the deep region (see Fig[7||. 



C. Third numerical experiment 

In the third numerical experiment [referred to as (iii)], 
the following function Z{(^, t) was used: 

Z(C, t) = C + A{\t)^e-^' In [1 + C exp(iC)] , (39) 

with A = tt/IOOO, A 0.25, C 0.95. As previously, 
p{d,Q) = 0, x(i?,0) = 0, and a(0) = gp. The results of 
this simulation are presented in Figs l8ll(7l 
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FIG. 10: (iii) t = 11.40 
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IV. SUMMARY 

In this paper exact nonlinear equations of motion for 
potential 2D ideal incompressible flows with a free surface 
in a uniform gravitational field have been obtained, tak- 
ing into account a time-dependent curved bottom with a 



profile prescribed by an analytical function. The math- 
ematical theory of functions of a complex variable and 
conformal mappings was extensively used in the deriva- 
tion. Despite a complicated structure of the equations, 
an efficient numerical code has been developed for space- 
periodic solutions, which gives very accurate results. 
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